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Abstract 

Using flow equations, equilibrium and non-equilibrium dynamics of a two-level system are in- 
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1^ I vestigated, which couples via non-commuting components to two independent oscillator baths. In 



equilibrium the two-level energy splitting is protected when the TLS is coupled symmetrically to 
both bath. A critical asymmetry angle separates the localized from the delocalized phase. On the 



'^ i other hand, real-time decoherence of a non-equilibrium initial state is for a generic initial state 

O ' faster for a coupling to two baths than for a single bath. 
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I. INTRODUCTION 

Under the notations of frustration of decoherence or quantum frustration effects are sub- 
sumed which are ascribed to the competition and mutual cancellation of two environments, 
which couple to non-commuting observables of a central system. The notion was coined in^ 
and the effect has since then been studied in a variety of systems, like a two-level system 
(TLS) coupled to two oscillator bathi"- or to two spin-baths^, a harmonic oscillator coupled 
to two oscillator batb^"- in spin-lattices^ or Josephson networks^. Most notably it was pro- 
posed as cooling mechanisnt^. The relation to Kondo physics was already pointed out ini. 
Certain phenomena occuring in the two channel Kondo model or in the Bose-Fermi-Kondo 
model can actually be interpreted in terms of quantum frustration^iii^. 

In the model originally studied inii^, a TLS with energy gap A couples linearly with 
its two transversal components to two independent baths. It will be called 2BTLS in the 
following. The strength of the ohmic coupling is measured by two quantities 73 and 72 
(assuming a magnetic field in x direction, bath 1 couples to the z-component and bath 2 
the |/-component). One remarkable result of Ref.- were the renormalization group (RNG) 
equations 



H'' 



-27^7?' - 7^ 'A= 



dl 

^ = -27™7f-7fA'^ 

f^(l-7$-'-7f')'». (1) 

where dl = —dlnuc is the differential of the flow parameter and h = A/uc, where Uc is the 
cutoff frequency of the bath modes. If either 73 or 72 is zero, the RNG equations of the 
single bath spin-Boson modeli^^i^ are recovered which predict a Kosterlitz Thouless phase 
transition for 7 = 1. For 73 = 72 ^^e renormalization flow is different: h scales always to 
infinity, i. e. a phase transition never occurs, not even for arbitrary strong coupling. This is 
by now one of the most striking signature of quantum frustration. 

However the question whether for large couplings the delocalized phase at symmetric 
coupling is stable against asymmetries remains unanswered by the above RNG equations. 
They do not yield any estimate for the renormalized energy gap A^, respectively Kondo 
temperature in the delocalized phase. 



The body of publications, mentioned above focuses on thermal equilibrium. But the 
question whether or not quantum coherence of a non-equilibrium initial state is protected 
by quantum frustration is crucial for possible applications. Time evolution of a spin in non- 
equilibrium can be more complicated than exponential decay predicted by Bloch equations^^. 
In particular an initially decoupled central system might on a very short time scale, called 
quantum Zeno-time, incur initial slips. This happens for instance to the dissipative harmonic 
oscillator—. In this case short times decoherence is indeed enhanced by a second bath and 
only later effects of quantum frustration occur-i^ . 

We address the above questions for the 2BTLS using the method of Hamiltonian flow 
equations. Flow equations were introduced in the early nineties by Glazek and Wilson^^ and 
about the same time by Wegner— . The method rests upon a continuous diagonalization of 
the Hamiltonian, details can be found in^. It was applied to the single bath spin-Boson 
model in^""— . In particular it proved to yield good results for the renormalized energy gap 

In this work a generalization of equations ([1]) is derived analytically, which embraces any 
kind of coupling to two baths. Numerically A^ is calculated as a function of an asymmetry 
angle, called 9, which varies from zero (single bath) to -7r/4 (completely symmetric). Whereas 
for weak coupling there is little dependence on the asymmetry angle, as the coupling becomes 
stronger the dependence on the asymmetry becomes more and more important. A symmetric 
coupling protects the gap and prevents the KT-phase transition. Identifying the critical 
angle allows us to plot a phase diagram in the 73 -72 plane, where the localized and the 
delocalized phases are separated by a critical line. 

Using techniques developed recently2^>^ we address the question whether decoherence 
of a non-equilibrium initial state is protected by a second bath. The answer to this can 
not be given without a careful distinction about what is meant by quantum decoherence. 
In a folkloristic definition decoherence is the decay of the off-diagonal elements in some 
pointer basis and relaxation the decay of the diagonal elements. For a two-level systems 
both processes are obviously not independent and it is therefore not easy to distinguish 
them. 

For symmetric coupling we find that the moduli of off-diagonal elements in the eigenbasis 
of the spin operator in x-direction incur initial slips and subsequent oscillations on a time 
scale of the cutoff-frequency ojc- These initial slips on the time scale of Uc are absent for a 



single baths, however the subsequent decay is oscillatory also in this case. The expectation 
value of the spin operator in x direction behaves quite differently. Here the decay is initially 
faster for a single bath but slows down on the time scale of the Rabi frequency A~^. On the 
other hand for symmetric coupling the decay is initially slow but increases later to reach an 
equilibrium value, which is smaller than for a single bath. 

In the first two sections of the manuscript we set up the model derive the flow equa- 
tions and calculate equilibrium quantities. In section IIVI the non-equilibrium djTiamics is 
considered. 

II. FLOW EQUATIONS FOR THE 2BTLS 

The Hamiltonian of the 2BTLS is given by 
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Fo = -A5i + X^5^a;fai,,a„,, (2) 



Hi'^ = Ss®Yl ^a!! («i,^ + «U) + ^^2 ® 5^ Ag (a2„ 
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k k 



where Si are spin |-matrices and an,k are bosonic annihilation operators [a„_fc, a]^ ^,,] = Skk'^nm 
and [an,fc,flm,fc'] = Wnk^'^mk'] ~ ^- ^^ ^^^^ ^^'^ ^^^ "^0 = \^2- The sum runs over the TV 
bath modes, where A^ is assumed a large number such that the spectral functions 

jf\.) = Y.{X^}f6{.-ur^) (3) 

I 

of both baths are smooth functions. They obey an Ohmic power law for small frequencies 
J3 '{(jj) = 273 'uj and J2 (w) = 272 UJ and are regularized by a cutoff Wc ^ A. For simphcity 
we assume here and in the following the cut-off and the number of bath modes to be the 
same for both baths. 

The Hamiltonian is approximately diagonalized by a unitary transformationi^ii^ which 
depends continuously on a flow parameter B. Any one-parameter family of unitarily equiv- 
alent Hamiltonians obeys the equation 

j^H^'\B) = [^^'\BlH^'\B)] (4) 



with a properly chosen anti-Hermitian operator r]^'^\ If t]^^) is chosen as the commutator 



V 



(0) 



r(0) 



[Ho(B), Hj (B)] it can be readily shown, see e. g.— that if Hq is non-degenerate 
in the limit B — )■ oo, tr H(B)Hj — > and thus the Hamiltonian becomes diagonal. The 
commutator on the right hand side of equation (^ generates interaction terms not present 
in if (°^ They are formally included in a more general Hamiltonian H = H^'^^ +H^^^ and in a 
new generator rj = [Hq, H]. The equations are closed by neglecting normal ordered products 
of more than two creation or annihilation operators. 

In order to write the interacting part of the form invariant Hamiltonian if in a compact 
form it is useful to arrange the creation and annihilation operators in a AN vector A = 



[ai' ,[ai ')%a2' ,{aW), where a^ = (a„,,i, . . . , a„^Ar), (al)^ = {al^, . . . ,alj^), n = 1,2. 



t^^ 
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It turns out useful as well to introduce coupling constants A^ [. = /vg ^ 



An 7. i A 



them in a AN vector A 



A 



Moreover S± = {S3±tS2)/2. Then 



(1)T ^-.(l)r ^{2)T ^{2)T^ , ^{n)T ^ , („) 



, A_,_ , A . 



, where A ± '' =(A^''/i, . . . , X^^^)^ ^ 
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Hi = S+® K^A + S. ^A^A + Si® : A^TA 



(5) 



The symbol : ab : denotes normal ordering with respect to a thermal expectation value. The 
4N X 4N matrix T has the following block structure 

■Sll til ■§12 ti2 \ 
tu -Sll ^12 Si2 



T 



U 



ti 



^ii 



^ij'i '^ij ^ ^ 



■^12 ^12 '^22 ^22 
V ^12 -512 ^22 S22 / 

Note the invariance of T under the unitary automorphism T — t- S~^TSa;, where S^ 

1 



(6) 



I2 ® cr^ 
reads 



Itv and a^ 



1 



is a Pauli matrix. Likewise we define S^. The generator 



ri = S+®A^{A-n)A-S^®A\A-n)A+Si®: A\n,T]A 



,W^ 



(7) 



where fi = diag (a;*^^\ — a;(-^\cj^^\ — w^^^) and w^"^ = diag(a;| , ...,cj]y'''), n = 1,2. In for- 
mer treatments of the Spin-Boson model with a single bath^^^^ within the flow-equation 
approach, a formally simpler generator was used instead of the canonical one 77 = [Ho,H]. 
This reduced the number of differential equations to be solved. The different generators 
were contrasted in Ref.— . 



In general there seem to exist by now no other guidehne to improve the canonical gen- 
erator than educated guess or physical intuition. Thus, for the present problem we stick to 
the canonical one. 

The commutator [r], H] is calculated straightforwardly and a set of non-linear coupled 
ODE's is obtained for the tunnelling matrix element A, the couplings A and for the matrix 
elements of T. They read 

^ ^ - -A (A - fi) coth -J— ^ A , 



dB 2 ' ' V 2 ^ 

^ = -(A-fifA + {T(A-fi) + [fi,T]}coth(M^A, 

^^ = -[Q, [fi, T]] - ^A (A - Q) k" - is,.A (A - ^) A^E, . (8) 

The equations (IHl) form a set of 1 + AN + 2N{AN + 1) first order non-linear differential 
equations which must be solved numerically. Before we do so, we show how they reduce to 
the RNG equations ([2]) for an ohmic bath in the low frequency limit. We limit ourselves to 
zero temperature. The differential equations for entries of T are of the type 

^ = cof{B)+giB), uER (9) 

which can be solved exactly 

f{B) = /(O)e'^^ + / rffi'e"(^-^')(7(fi') . (10) 

Jo 

This might be plugged into the flow equation for A. It suffices to evaluate these equations 
for small frequencies. Using the definitions of the spectral functions (|3]) and 



E >^^kB)Xp^,{B')6{u - u;P) = 2^^f^fu , Vj, j' G {2, 3} (11) 



k 

an integro-differential equation for the coupling constants is acquired 
^^^^ = -2A{bM^\b) - 2 /''d5V7r(i?)7r (B'H 



2A{Bf^t\B) - 2 / dB'^J^t\B)ltHB')ul 
Jo 

■^0 m=l V ^ 

x{2A{B')^J^,t\B')-u^^J^t\B')x^ . (12) 
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The corresponding equation for 72 is obtained from Eq. (TT2|) by interchanging the indices 2 
and 3 everywhere. This equation allows for a perturbative expansion in h = A/uc- Keeping 
only the highest order term in the integral Eq. flT2|) reduces to 

X /' dxa;e-^--'(^-^') V ^J ir\B)lt"\B') . (13) 

In the limit Uc ^ 00 the B' integration becomes (5-like for almost all x G [0, 1] and we arrive 

at 

dlt\B) 



dB 



-2ulh\B)^t\B) - ^lt\BH E it'^B) (14) 



m=l 



and likewise for 72 (-B). To make contact with the RNG equations, we use the relation^^ 

1 



We = ^^ = e-^ (15) 



and obtain 



,W 



^=/z^7f-27f5:7r^ 

m=l 



d/ 



/^'7f -27^ E^2'"^ r. = l,2. (16) 



m=l 

One derives straightforwardly from equation ([8]) 

rf/i(/) ^ ^ 



EE^rh- (17) 



\ j=2,3m=l 

Equations ( TT6l) and ( IT71) correspond to the one-loop perturbative renormalization group 
equations for arbitrary couplings 73 and 73 , n = 1,2. We do not analyze them further 
here, but only mention that the result of Novais et al.- stated in Eq. ([1]) is obtained by 
setting 72 and 73 to zero. However it must be pointed out that the same equations are 
obtained for 73 and 72 5 if Ts and 73 are set to zero, i. e. in the absence of the second 
bath. 

An adaptive step-size fourth order Runge-Kutta algorithm has proved to be a reliable 
solver of the flow equations ([8]). Most entries of A and of T become exponentially small for 
large flow parameter and the Hamiltonian becomes diagonal 

H{oo) = -ArSi + ^ : A^\Q\A : +H,,, (18) 
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FIG. 1. Left: Plot of the renormalized tunneling matrix element A^ as a function of the angle 9 
defined in the main text, the total coupling strength is 7tot = 0.1 (crosses, online yellow), 7tot = 0.3 
(empty circles, online orange), 7tot = 0.5 (filled boxes, online red), 7tot = 0.8 (empty boxes, online 
dark red) and 7tot = 1 (filled circles, full black line). The cutoff frequency is ujc = lOA. The 
number of bath modes is A^ = 1000. 
Right: the same for the equilibrium expectation value (^i). Uc = lOA, N = 400. 



with a finite renormalized tunnelling matrix element A^ = A(oo). Not all entries of A and of 
T decay exponentially for large B. From the flow equations ([8]) it is seen that the coupling 
matrix elements X^[. for frequencies close to the renormalized tunnelling matrix element 
decay most slowly. On the other hand the diagonal entries of T do not decay at all, leading 
to an effective coupling of the bath modes to Si in the renormalized Hamiltonian 

2 

Hrcs = 'S'l ® ^ ^ Snn,kki^)(^li,k^n,k ■ (19) 

n=l k 

Although this term - being diagonal - causes no additional difficulties, for practical purposes 
it can be neglected, since the residual matrix elements Su^kk, S22,kk are usually much smaller 
than the mean level spacing of the bath modes. 

In figure [1] the renormalized energy gap of the two-level system is plotted for a fixed 
overall coupling 7tot = 73 +72 ^is a function of the relative angle 9 = arctan I -^72 773 J 
which varies from zero (single bath) to 7r/4 (equal coupling strength). Whereas for small 
overall coupling the renormalized energy gap A^ is almost independent of 6, for increasing 
coupling strength the gap is protected by a symmetric coupling. Finally for 7tot = 1 the 
energy gap renormalizes to zero for 6* = but remains finite for symmetric coupling. 
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FIG. 2. Phase diagramm in the 7tot ~ & plane. The line indicates the critical asymmetry angle, 
which separates the localized from the delocalized phase. The critical angle was determined for 
N = 800 bath modes. 



If 7tot is increased even further the energy gap A (5) crosses zero for some large value 
of B and decays afterwards very slowly in an oscillatory fashion to zero. This happens 
for angles smaller than some critical angle, indicating the onset of the strong coupling 
regime, respectively of the KT phase transition. It is expected that the flow equations, 
being generically perturbative, become less exact for stronger coupling. However for 6* = 
the critical value 7 = 1 was obtained analytically and with good precision numerically^ . 
Therefore it is well justified to assume that the flow equations yield a good estimates for the 
critical 7tot for 7^ as well. 

In figure [2] the critical line is plotted in the 7tot ^ & plane, which separates the localized 
from the delocalized phase. It is seen that it crosses the x-axis at some value smaller 
than one. This offset is due to the finite number of bath modes and of the finite cutoff 
frequency. This can be improved systematically by increasing the number of bath modes 
and simultaneously increasing the endpoint of the flow i^max- For values of 7tot larger a than 
some value 7tot ~ 2.5 the flow becomes unstable. 



III. EQUILIBRIUM EXPECTATION VALUES 

In oder to calculate equilibrium expectation values with respect to the transformed Hamil- 
tonian H{oo) the corresponding operators have to transform as well. Complex 4A^-vectors 
X and C 1 ^^^ introduced and the spin operators are expanded as 

^1 = hoSo + hiSi + S+® X^A + S- ® A^x 

s+ = h+s+ + h^S- + 5o ® c U +Si® cIa 

S^ = h*^S+ + hlS^ + So ® A^Co + Si ® i^C 1 (20) 

The flow equations for ho, hi, h± and for x ? C o i ^^^ obtained by calculating the commutator 
[77,5*4]. The equations are closed by neglecting all normal ordered operator products with 
two or more annihilation or creation operators. They are stated in App. [Bl The equilibrium 
density matrix with respect to the renormalized free Hamiltonian flTSl) is just 

Peq = (so + tanh (M^ Si^ ® pS ® pi? . (21) 

Here piq = Y[k^^P(~l^^k 'K.k'^rL,k)/ is the thermal density matrix of the two free envi- 
ronments. Thus, once the equations are numerically solved, an arbitrary equilibrium ex- 
pectation value of the spin operators is readily calculated. As an example we consider the 
one-sided Fourier transform 

XU^) = -W ^e^^'{[S3{0),Ssm (22) 

Jo ^TT 

of the correlator ([5*3(0), 53 (t)]) which was investigated in^. At zero temperature, its imagi- 
nary part x'zz is given by 

x'Lico)^{h+ + h.fSiu-Ar)+[Co + Ciya + ^.)H^-M){Co + Ci) ■ (23) 
As a second example, we consider the equilibrium expectation value 

(.,.|4.„k(M). (,4) 

It is plotted in the bottom picture of figure [T] for zero temperature and for different angles 
6 as defined before. Since the calculation is numerically more expensive than that of the 
energy gap, the number of bath modes is A^ = 400. For small and intermediate coupling 
it behaves qualitatively similar to the renormalized two-level energy gap A^. For strong 

10 



coupling 7tot ~ 1 it is seen that {Si) does not scale to zero for ^ = as expected, indicating 
that the flow equations lose accuracy in the strong coupling regime. 

Before we discuss the numerical results for the equilibrium correlation functions an ex- 
planatory remark is in order. A careful treatment of equilibrium correlation functions within 
the flow-equation approach requires high sophistication. For frequencies close to the renor- 
malized tunnel matrix element A^ the flow converges only very slowly with i?max, the end- 
point of the numerical integration of the flow. Since the endpoint of the integration is itself 
limited by the density of the bath modes an accurate resolution would require an out of 
scale number of bath modes. As a consequence of this numerical limitation the equilibrium 
correlation functions have a two-peak structure: one broad maximum at a value smaller 
than Ar and a second sharp peak right at A^, which is clearly unphysical. 

The problem can be overcome by employing constants of motion under the flow. This 
was done in2^ for the one-bath spin Boson model. The result is a smooth curve with a single 
peak. But such constants of motion under the flow are not always easy to identify. 

We refrain from this procedure and show the curves for x'Li^) obtained by fitting the 
numerical data with smoothing splines using an extremely high fidelity factor (of order 10^) 
everywhere but around A^, where it is quartically suppressed. 

In figure [3] the correlation function x'Li^) i^ plotted for equal coupling strength to both 
baths and with an overall coupling strength 7tot varying between 0.1 and one. The curve 
corresponds to Fig. 4 in reference^ and is qualitatively similar. As the coupling strength 
increases the resonance peak becomes smaller and smaller but never disappears. The max- 
imum of the resonance peak is systematically below A^.. This is a difference to Fig. 4 in 
reference^ where the maximum seems to be always right at the renormalized tunnel matrix 
element. 

In figure |3]the correlation function x'Li^) is plotted for fixed overall coupling strength 7tot 
and for different angles 6. The resonance peak in the symmetric case {6 = it /A) is largely 
enhanced as compared to the highly asymmetric case [6 = O.Itt). However, the reason 
for this is rather trivial. In the highly asymmetric case the coupling to the ^-component 
is largest, whereas there is no coupling to the y-component. In the symmetric case the 
coupling to the z-component is reduced, which is reflected by the enhanced resonance peak 
of x'zz- However the coupling to the y-component is larger, which yields a reduced resonance 
peak of Xyy i^ot shown here). If we write x'Li^^^) ^^ ^ function of the relative angle 9, 
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FIG. 3. Plot of the transverse susceptibility Xzz(w)/w in 2;~direction for symmetric coupling 73 
= 72 and for ten different values of 7tot = v2 • O.ln, 1 < n < 10, from top to bottom (online 
color: from dark-colored to light-colored). The number of bath modes is 400, A/wc = 1/10. 

then the obvious relation x^'^(c<;,^) = y^yyiuj.-n jl — 6) holds. Thus an enhancement of the 
resonance peak in z-direction comes necessarily with a decrease in y-direction and vice 
versa. Indeed in Fig. [3] the resonance peak of xL is biggest for 6 = O.Stt in spite of the 
asymmetric coupling (for even higher 6 it increases more and more). Note however that the 
location of the maximum of the peak is maximal in the symmetric case. 



IV. THERMALIZATION AND DECOHERENCE 



In thermal equilibrium the mutual energy transfer from the system to the environment 
and vice versa is zero, warranted by fluctuation dissipation theorems. However in the process 
of thermalization the net energy transfer of the system to the environment is positive. 

Assuming an decoupled initial state, which is fully polarized in some direction perpen- 
dicular to the X axis (we may assuume < 9' < 7r/4) 

Pinit = (5-0 + COS ^'^3 + Sin ^'^2) ® p« ® pg) , (25) 

thermalization is characterized by the time evolution of the expectation value of the system's 
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FIG. 4. Plot of the transverse susceptibility \'!^.^{(jj) / uj in z-direction for three different angles 9 
= O.lvr (full line), 9 = 7r/4 (dashed line) and 9 = O.Svr (dotted line) for overall coupling strength 
7tot = 0.3. The number of bath modes is A^ = 400, tJc/A =10. 

energy {Hs(t)) = —A{Si{t)) . This quantity is expected to approach its equilibrium value 
on a certain time scale, the so called relaxation time, which is usually denoted Ti. 

Decoherence is the creation of entanglement of the system with the environment. It is 
measured by the decay of the off-diagonal elements of the reduced density matrix of the 
spin in the Si basis, i. e. by the expectation values {S±). A basis independent measure for 
decoherence is the purity P(t) = P||(t) +V±{t), where P||(t) = 2 ^^^^(^^(t))^ and 'P±{t) 
= 2 ^,^^2('S'n(^))^- Decay of decoherence takes place on a time scale T2, called decoherence 
time^, we asociate it with V±{t). Both decoherence time and relaxation time enter in the 
definition of purity. We call the two quantities V±{t) and V\\{t) transverse respectively 
parallel purity. For the initial state ([25]) V±{0) = P||(0) = 1/2. 

Assuming a decoupled initial state as in Eq. ( I25l) first order differential equations for the 
spin expectation values are straightforwardly derived in second order perturbation theory 



Jt^'^^ 



Tf\t) + Tll\t)){Si)-F{t) 



d 
di 



.(1) 



{s,) = A),'>{t){Ss)-n'>ms2) 



Wi 
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|(53) = -Af(t)(52)-rf(t)(53), (26) 



with the time-dependent coefficients 



1(7*) (t) = dt' dcocos (A(t' - t)) cos (tuit - t')) Jt\co) coth {uj(3/2) 





t 00 



A([")(t) = A- /"rft' f dcosm{A{t'-t))cos{u{t-t'))Ji^\co)coth{u(3/2) 



t 00 

fW^/*'/<^s,„M.-.'))sm(A(.'-0)(4=>M + 4"H) . (27) 



In the Markov approximation these coefficients become time independent r„ (t) = r„ 
= (7r/2)J^^(A)coth(/3A/2), F{t) = F = (7r/2)(jf (A) + 4'\A)) and At\t) = A^^ = 
A — A J^" dco coth(/3aj/2) Jn {(^)/{uj'^ — A^). Note that for an ohmic bath and at zero 
temperature A„ has a logarithmic singularity in the cutoff frequency coc- 

From equations (!26l) the phenomenological Bloch equations are obtained which predict 
an exponential decay of decoherence and of relaxation. Their solutions are 



{Sr{t)) = {{S,m - {SiU)e-^'^ ^'^ ^' + {S,) 



(rr^+r«)t 



cq 



Wn{tj) - X - \ \ - X ^ ^ 

where A± are the roots of the characteristic polynomial 

xiu) = u^- rujivf + r(^)) + Af^AW - rf r(^) (29) 

Thus decoherence and relaxation time are given by Ti = l/(r2 + Fg ) and T2= 2Ti. In 
second order perturbation theory the friction coefficients of the two baths add up. No 
frustration occurs. 

In the Markov approximation Bloch equations hold beyond perturbation theory with 
relaxation and decoherence times depending in a more complicated non-perturbative way 
on the coupling strength 73 and 72 . Corrections were calculated in Ref.-. In the regime 
where the Bloch equations f l2B]) hold, the quantum regression theorem can be invoked and 
the dynamics of the expectation values is governed by the equilibrium correlation functions. 

Yet at low temperature and on the time scale of the inverse cutoff frequency, Bloch 
equations do not hold. The coefficients in Eq. fl26|) become time dependent and the simple 
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exponential behavior (128|) breaks down. This is seen most directly in a Taylor expansion of 
the time-evolution operator U{t) = 1 - iHt - t^H^/2 + 0{f). For the initial state (ES]) it 



,(1) 



,(2) 



predicts a quadratic behavior of (5*1) = t^/2r^ + 0{t^), where r ^ ^ ujcsj 'y^''' + "^2' i^ the 
inverse quantum Zeno time. For the transverse purity one obtains 



2P^{t) = 1 - — {sin\e) cos2(0') + sin\e) cos2(^')} + o{f) 



(30) 



The quadratic time dependence vanishes iff ^ = ^' = 0. This indicates that initially, for 
short times, a symmetric coupling accelerates decay of coherence. 

In order to monitor the time evolution of the expectation values in the transient regime 
on a time scale of order of the quantum Zeno time, methods of non-equilibrium real time 
thermodynamics must be employed. Real time quantum evolution is addressed within the 
flow equation approacb^ii^^ by applying subsequently the unitary transformation Ub{Bi, B2) 
generated by ri{B) and the time evolution operator Ut.oo{ti^t2) = e"*-'^'^"")^*^"*^) on the oper- 
ator of interest according to the diagramm: 

Ub{0,oo) 
O(B = 0,t = 0) , C(oo,0) 



0(0, t) 



Ub{0, -00) 



Ut,ooiO,t) 



0{oo,t) 



Since time evolution is simple for 5 = 00 the observables are first transformed into the 
B = 00 basis evolve in time and are then transformed back. At time t the Heisenberg 
operators have been propagated by the diagonalized Hamiltonian fITH]) . This yields new time 
dependent expansion coefficients h±{t)= e^*^''*/i±(oo), xit) = e**^^''"^)* x(oo) and Cn(^) = 
e*^*Cn(oo), n = 0, 1. The coefficients ho = ho and hi = hi remain constant under time 
evolution. These coefficients are numerically transformed back, yielding an approximate 
solution of the Heisenberg equation for the spin operators. The expectation value with 
respect to the density matrix ( l25l) are 



(Slit)) = h,{t)/2 

(52(t))=Im[(/i+(t)-r(t))e^^72 
{S,{t))=Re\{}i+{t) + h*_{t))e'''/2 



(31) 
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FIG. 5. Right: Time evolution of the transversal purity for an initial state characterized by the 

angle 6' = Q for different triplets (0,7tot,'^c/^)- These are (7r/4, 0.3, 10) (online black, crosses), 

(0,0.3,10) (online yellow, asterisks), (0,0.15,10) (online light red, boxes), (7r/4, 0.6, 10) (online 

dark purple, dots), (7r/4, 0.3, 20) (online lighter purple, triangles), (0,0.3,20) (online darker red, 

triangles). 

Left: The same for short times on a logscale. 

The calculation is numerically delicate^"^. In order to perform the backward integration 
the forward flow of the Hamiltonian must be stored. This is a sizable amount of data of order 
of one terabyte. The read-in and the read-out slow down the routine. We thus performed 
the calculation of V_[_ with 250 bath modes, respectively of {Sx{t)) with 100 bath modes. 

In figure [5] the transverse purity is plotted for different values of 73 and 72 for an 
initial state characterized by the angle 6' = 0. It is seen that for short times of order of 
u~^ the transverse purity decays faster for a symmetric coupling than for a single bath, as 
predicted by equation (1301) . The decay occurs in an oscillatory fashion for both a single bath 
and for symmetric ccoupling. Although the dissipative two-level system has been studied 
extensively^^ to our best knowledge this oscillatory purity revival was not reported before. 
By now we do not have a satisfactory physical explanation for it. For symmetric coupling 
the oscillations decrease rapidly in less than one period of the Rabi oscillations. As can be 
seen from left picture of Fig. [5] the frequency seems to scale with Uc and the amplitudes with 
7tot- For a single bath the oscillations are much slower and decay less rapidly. 

The dependence on the initial state is considered in Fig. [61 The transverse purity for the 
initial state characterized by 6*' = and for the initial state 6' = n/4 is plotted. Whereas for 
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FIG. 6. Time evolution of the transversal purity for an initial state characterized by the angles 
9' = 0, 7r/4 for symmetric coupling (asterisks, online black) and for a single bath, 0' = (dotted, 
online red) and 9' = 7r/4 (boxed, online orange). The other parameters are 7tot = 0.3 and uJc/ ^ 
= 10. 

symmetric coupling there is no visible difference, for a single bath the initial decay is much 
faster for 6' = tt/A than for 6' = 0, see Eq. (|30|). 

The time evolution of {Sx{t)) is plotted in Fig. [7] for symmetric coupling (9 = 7r/4) and 
for a single bath [6 = vr/4) for a moderate overall coupling strength 7tot = 0.3. Here the 
expectation value indeed decays initially faster for a single bath than for symmetric cou- 
pling. However on a time scale of the Rabi-oscillations the decay grows faster for symmetric 
coupling to reach an equilibrium value, which is smaller than for a single bath in accordance 
with Fig. [H 

V. SUMMARY &i DISCUSSION 

While the calculation of equilibrium correlation functions is somewhat cumbersome within 
the flow equation approach, the method turns out to be a useful numerical tool in non- 
equilibrium physics. We were able to monitor purity decay on the time scale of the quantum 
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FIG. 7. Time evolution of the expectation value {Sx{t)) for the initial state (I25p for symmetric 
coupling (online blue, stars) and for a single bath (online red, crosses), 7tot = 0.3, uJc/A.= 10. 



Zeno time as well as on the time scale of the inverse Rabi frequency. 

When one speaks about coherence of a two-level system one has carefully to distinguish 
between the decay of the off-diagonal elements and of the diagonal elements. It is char- 
acteristic for a small size Hilbert-space that both are not independent and the distinction 
between decoherence and dissipation is fuzzy. 

In our analysis frustration effects of two independent oscillator bath could be identified 
in the renormalized energy gap A,., in the ground state expectation value of 5*1 and in the 
ground state energy shift. These quantities are protected by a symmetric coupling. In 
particular the protection of (^i) can rightly be called protection of decoherence since it 
contributes to a high equilibrium purity of the spin. 

In non-equilibrium relaxation, i. e. the decay of (Hs) oc — (^i), is protected by a sym- 
metric coupling on a time scale of the quantum Zeno time. However, the decay of the 
off-diagonal matrix elements of the reduced density matrix, corresponding to (5*2), (S^) and 
to the transverse purity is systematically faster for a symmetric coupling. 

The decay of both (Hs) and of the transverse purity occurs in an oscillatory fashion. The 
physical reason behind these oscillations is unclear. 
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The dependence of the renormahzed energy gap A,, on an asymmetry angle is a generic 
non perturbative effect. The flow equations ([8]) might be truncated by setting all second 
order terms, i. e. the matrix entries of T (Eq. (jSD), to zero. The truncated flow equations can 
be analyzed analytically, see App. \M The outcome is ln(Ar/A) oc — 7tot/(l — Ttot), similar 
to the old result by Silbey and Harris2£ which features no dependence on the asymmetry 
angle. Our analysis affirms that the delocalized phase for couplings 1 < 7tot < oo is stable 
against small asymmetries. 

The perturbative RNG equations (1161) and (11 7p obtained from the flow equations are 
completely sjTumetric in the four coupling constants 72 ? 73 ? 1^ = 1)2. Setting any two of 

them to zero yields the RNG equations of Ref.-, with the implication of a delocalized phase 

(2) (2) 
for 7tot -> 00. Setting for instance 72 =73 =0, this implys that also a symmetric coupling 

of the spin with its y and z components to a single bath can protect the delocalized phase. 

This question requires further investigation. 
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Appendix A: Linearized Flow equations for two baths 

We consider the linearized version of the flow equations. In the linearized version of the 
flow equations the flow of T can be neglected. 

For ohmic spectral functions J^'^'{u) = 2'yl'^'u9{uc — u) , -i = 2, 3, n = 1, 2 immediately the 
first order RG equations 

^ = -A%f) z = 2,3, n = l,2 (A2) 
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are obtained. Introducing the auxiliary densities 

Gi")(a;) = J2 (^£i)' ^ (^ - 4"0 ^ = 1' 2 (A3) 

k 

the renormahzation group equation for the tunnehng matrix element ([H]) can be written as 
dA{B) I'^v- /"j coth(/3a;/2) d ,.. ^. ,,^, 

n=l (T=± "^ ^ ^ 

Following the outlines oi^ a self consistency equation for A^ can be obtained. For zero 
temperature it reads 

l„^^VV/^fV.°) (A5) 

A ^-^ ^-^ \ 4A,. A,, - o-w ^ ^ 

n=l a=± ;^ ' ' 

For A^^i = A^^i = 0, di^ = G'l^ = 4^^ and Cf = gL^^ = jf ^ and for an ohmic bath the 



renormalized matrix element becomes 



.i'h.i'^ 



A, = A(|)^^^r^. (A6) 

This is a straightforward extension of the old result by Silbey and Harris22. In the linear 
approximation of the flow equations there is no angle dependence of A^,. The full flow 
equations must be employed. 

Appendix B: Flow equations for the spin operators 

The flow equations for the expansion coefficients of the spin-operators are obtained from 
the commutators [77,5*1] and [?7, Si]. They read: 

dB 2 ^ ^ "^ 

^ = -2 A (A-fi)coth^^jx 

''^^^^-Uc:cothrMVckU^,(A_^)A 



dB 2 V V 2 



dB 2 V V 2 

^ = h{A - m + [n,T] coth (^ ) X 

dUB) 



dB 

dCiiB) 1 



[fi,T]S,Ci 

^^ 2 (^- + ^+^^) (^ " ^^^ + [^'^]^^^Co (Bl) 
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These differential equations are the same for the forward flow and for the backward flow. 
However the initial conditions are different. For the forward flow the initial conditions are 
hi = h^ = 1 and all other components are zero. Since the differential equations are linear 
in the expansion coefficients the imaginary parts of h±, x and Cq ^ remain zero throughout 
the flow. 

Due to the time evolution the imaginary parts acquire a non-trivial backward flow. The 
initial conditions are now Rex (t, 0) = cos((r2 — Ar)t) Rex (oo), Imx (t, 0) = sin((r2 — Aj.)t) 
Rex(oo), ReCo,i(^!0) = cos{Qt) ReCo,i(oo) , Ini^o,i(^)0) = ReCo,i(oo) sin(r2t), Reh±(t,0) 
= cos(Art) h±{oo) and Im/i±(t, 0) = ±sin(Aji) Re/i±(oo). The flow of the imaginary parts 
Imx decouples from that of the real parts and of /iq and of Hi. Thus it needs not to be 
considered. 
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